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Abstract 



We calculate microscopically the properties of the dilute 3 He component in 
a 3 He- 4 He- mixture. These depend on both, the dominant interaction between 
the impurity atom and the background, and the Fermi liquid contribution due 
to the interaction between the constituents of the 3 He component. We first 
calculate the dynamic structure function of a 3 He impurity atom moving in 
4 He. From that we obtain the excitation spectrum and the momentum de- 
pendent effective mass. The pole strength of this excitation mode is strongly 
reduced from the free particle value in agreement with experiments; part of 
the strength is distributed over high frequency excitations. Above k > 1.7A~ 4 
the motion of the impurity is damped due to the decay into a roton and a low 
energy impurity mode. Next we determine the Fermi-Liquid interaction be- 
tween 3 He atoms and calculate the pressure- and concentration dependence 
of the effective mass, magnetic susceptibility, and the 3 He- 3 He scattering 
phase shifts. The calculations are based on a dynamic theory that uses, as 
input, effective interactions provided by the Fermi hypernetted-chain theory. 
The relationship between both theories is discussed. Our theoretical effec- 
tive masses agree well with recent measurements by Yorozu et al. (Phys. 
Rev. B 48, 9660 (1993)) as well as those by R. Simons and R. M. Mueller 
(Czekoslowak Journal of Physics Suppl. 46, 201 (1996)), but our analysis sug- 
gests a new extrapolation to the zero-concentration limit. With that effective 
mass we also find a good agreement with the measured Landau parameter 
F a - 



I. INTRODUCTION 

Many ground state properties of 3 He- 4 He mixtures, in particular the energetics of the 
system and the local structure, are today quite well understood both experimentally^ and 
theoretically from a microscopic point of view]! With "microscopic" we mean that one 
postulates no more knowledge than the empirical Hamiltonian 
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^ = -EE^+ E E niri a, -rf |) (1.1) 

a i=l a a/3 i,j 

that contains only a local two-body interaction (recent work uses most frequently the AziJ 
interaction) and the masses of the two species of particles. This paper builds on our pre- 
vious studies of ground state properties of 3 He - 4 He mixtures^ and properties of single 3 He 
impuritiesB in 4 He. These calculations have produced an overall accuracy of the energet- 
ics of the order of a few tenths of a percent for all experimentally accessible densities and 
concentrations. 

To specify our notation we use in the following Greek subscripts a, ft, . . . G {3, 4} to refer 
to the particle species (a 3 He or a 4 He particle), and Latin subscripts i, j, ... as in the to 
refer to the individual particles. The prime on the summation symbol in Eq. ( |1.1|) indicates 
that no two pairs (i,a), (j,ft) can be the same. The number of particles of each species is 
N a , and iV = N 3 + N± is the total number of particles in the system. In terms of the 3 He 
concentration x we have 



N 3 = xN, N 4 = (l-x)N (1.2) 

and the corresponding densities = N a /Q are inversely proportional to the volume fl 
occupied by the whole fluid. 

The overall most successful microscopic theory of the helium liquids is the Jastrow- 
Feenberg variational method&i which uses a variational ansatz for the ground state wave 
function of the form 

^o({r l (a) })=e^^ ) »$ ({rf)}) 

N a ,N N a ,N fj ,N 1 

U({r\ a) })= ^T, E' « (< *W;) + ^j£ E' « M7 W,,r fc ). (1.3) 

a/3 i,j ' q/?7 i,j,k 

Here $o({ r i }) is the Slater determinant of plane waves ensuring the antisymmetry of the 
Fermion component of the wave function. The functions u^ a ^\vi, Vj) and u <ya ^\v i , Tj, r^) are 
the pair and triplet correlations; the species superscripts determine the type of correlation. 
An essential part of the method is the optimization of the ground state correlations by the 
variational principle M 

SE ° 0, ., =0, (1.4) 



where 

is the variational energy expectation value. Details of the procedure, and the necessary 
working formulas have been discussed at length in Ref. 0. 
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II. DYNAMICS OF A SINGLE 3 HE ATOM 



This section briefly reviews our methodSS to calculate the dynamics of an impurity 
atom in the host liquid 4 He in its ground state with the wave function ^(ri, ...,r^). The 
Jastrow-Feenberg wave function for the ground state of a single impurity atom in 4 He is 
easily derived from the wave function fll.3|) by setting A^3 = 1 and the Slater determinant 
$o( r o) equal to unity, 



i N i N 

* (3, (ro,ri rjy) = exp - [$>( 34 )(r , r,) + - £ u( 344 )(r , r if r fc )] *(r lf r N ) . (2.1) 

j=l ' j,k=l 

The coordinate ro refers to the 3 He atom and ri . . . rjv to the 4 He atoms. Results for ground 
state calculations with this wave function have been reported elsewhere.B They are consistent 
with the low concentration limit of the calculations of Ref. ^| and in quantitative agreement 
with experiments. 

The dynamics of an impurity atom is determined by its response to a weak, external 
time-dependent perturbation £/ ext (r ; t). A natural generalization of the wave function ( [2.1| ) 
for a moving impurity atom is to allow for time- dependent correlations. The kinematic and 
dynamic correlations are separated by writing the wave function in the form 



$(t) = e -Ufcr+i*/» #( 3 )(r , r x , ...r N ; t)/yJ(W(t) \ W{t)) , (2.2) 

where E^ + i is the variational ground state energy of the N + 1 particle system, and 
\l/( 3 )(ro, i"i, ...rjy; t) contains the time-dependent correlations, 

1 N 

*®(r , n, ...r N ; t) = exp - [6u®{r ; t) + £ Su^\r , r,; t)] ¥ 3 \r , n, r N ) . (2.3) 

1 i=i 

The time-dependent components of the wave function are determined by an action principle, 
searching for a stationary value of the action integral 

C = Jdt($(t)\H^ + U cxt (r -,t)-ih-mt)) , (2.4) 

where is the Hamiltonian of the impurity-background system, obtained from Eq. ( |1 . 1| ) 
for N 3 = l. 

We make two assumptions in the evaluation of the action integral. First we require 
that the pair- and triplet correlation functions in the ground state are optimized. This 
is important because it eliminates all contributions to the action integral ( |2.4|) that are 
linear in the time- dependent correlation functions. Then we assume that the perturbation 
is weak which allows us to keep only the quadratic terms in 5u^(r ;t) and 8u^ u '(r , r$; t). 
The expression Q2.4p simplifies because the potential energy term commutes with the time- 
dependent part of the wave function and we are left with second-order terms originating 
from the kinetic energy and the time derivative only.00 

The variation of the action integral with respect to 5u^(vQ]t) and 5u^(r , r^; t) leads 
to one- and two-particle continuity equations, respectively,0'i 
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V • j (3) (r ;t) + pp [3 \r ;t) = ^f/ cxt (r ;f) 



(2.5) 



V • j (34) (ro,r i; t) + Vi • J (34) (r , r l5 1) + J^Vo, ri; = \u^{v Q - *)p (34) (r , n) . (2.6) 

We have kept only the time dependent parts of the full densities. The transition currents 
are defined in terms of the fluctuating one-particle density and pair correlation function: 



f\r ;t)- 



h 



2m 3 i 



V *p (s) (r„;f) -fd'n Su^(r , n; i)V p (34) (r , n)l , (2.7) 



j (34) (ro,r 1 ;t)=4)j (3) (ro;t)p (34) (ro,r 1 ) + 



n 



+ J d 3 r 2 



p( 344 )(r , ri ,r 2 )- 



2m 3 i 



p( 34 )(r ,r 1 )Vo^ 34 )(ro,r 1 ;0 



,(3) 



p( 34 )(r , ri )p (34) (ro,r 2 )lV 2( 5 M ( 34 )(r ,r 2 ;t) 



J (34) (ro,r i; t) = 



ft 



2m 4 i 



,(34), 



ro.rOVi^Hro.ri;*) . 



(2.8) 
(2.9) 



/From the calculation of the ground state properties of the impurity we need here the 
impurity-background pair and triplet distribution functions p( 34 )(r , r x ) and p( 344 )(r , ri, r 2 ), 
respectively. 

We assume that an infinitesimal external potential drives the system with a given fre- 
quency and wavelength, and the system responds with a density fluctuations of the same 
frequency and wavelength. This determines the linear response function 



(2.10) 



The inverse of the response function can be calculated from the Fourier transform of the 
continuity equations. The Fourier transform of the density fluctuations is defined as follows 



5p {3) (k, u) = J d 3 r dt e -iO*-*°-»*)SpW (r ; *) , 



(2.11) 



and analogously for the external potential. 

The poles of x^(k,u) determine the elementary excitations; their dispersion relation 
is obtained from the first continuity equation by setting U cxt (k, u) = 0. This leads to the 
implicit equation 



h 2 k 2 



with the self-energy 



^ 3 \k,u) 



h 2 



2m 3 



d 3 p 



2m 3 J (2tt)V 4 ) \hw-tW(k + p) -e( 4 )(p) 



(2.12) 



(2.13) 



Here S^ u \p) is the 3 He- 4 He structure function, t^(k) the kinetic energy of the impurity 
and e^(p) the background phonon-roton spectrum. The function (3^ (p) will be defined 
below. 
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The linear response function has then a simple form 

X {3) (^) ' 



HOU - (k) - S( 3 ) (k, U) + 17] hu + t(3) (jfc) + S( 3 ) (Jfe, w) + Z7? ' 



7] -> + 



(2.14) 



This is the density-density response function of the 3 He component in the dilute limit. The 
response function is related to the Green's function throughtJ 



d 3 q d{%u') 

(27T) 4 



G(q,u') G(k + q,oj' + Oj) 



where 



+ 



hu-t^{q)-^){q 1 u)+i7 1 hu-t^){q)-^){q 1 uj)-i7 1 ' 



(2-15) 



r]^0+ (2.16) 



and n q is the Fermi distribution. In the dilute limit, we have n q 



5 q and Eq. (ETTBj) 

reduces to Eq. (|2.14j) or, vice versa, the first part of Eq. ( |2.14j) can also be regarded as the 
single-particle Green's function. 

The imaginary part of the linear response function defines the dynamic structure function 



S(k, uj) 



%m [x i3) (k,u)] , 



7T 



(2.17) 



which can be measured in scattering experiments. 

The contributions of the elementary excitations can be separated out as a delta-function, 



S(k, uj) = Z(k)8(huj — Twjp) + S m (k, uj) . 



(2.18) 



Here Tiuj p is the solution of Eq. (|2.12|) and S m (k,u) is the contribution of the modes that 
make the self energy complex. This becomes possible when the denominator in Eq. (|2.13 ) 
becomes positive for some value of p, 



max 
p 



^-t (3) (k + p) -e^(p) 



(4), 



> . 



(2.19) 



If the condition (|2.19 ) is satisfied, then it is kinematically possible that the 3 He impurity 
loses energy by emitting a phonon-roton mode (p) while making a transition into a low- 
energy impurity mode t^ 3 )(k-(- p). xhe strength of the pole Z(k) can be evaluated from the 
derivative of the self energy, 



Z(k) 



d(huj) 



(2.20) 



The second continuity equation ( |2.6| ) can be written in the Fourier space and formulated 
as a linear integral equation for /3^(p), which is related to the fluctuating pair correlation 
function as discussed in the Appendix. The integral equation sums ladder diagrams up to 
infinite order, 
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hu 



k-p S {u) (p) r d 3 q 

SM(pj ~ J (2tt)3 p (4) \ hu _ t(3 )( k + q) 



,(4) 



(?) 



(2.21) 



with the kernel 

# kl .(p, q) = ^ (44) (^){[(^ (34) (|P - q|) + l) ^ (344) (P " q, "P, q) + S (34) (|p 
(k + p)-(k + q) fffl 



q|) 



hu — 



v 



^ (34) (|p-q|)^ (4) (p)- (2-22) 



Here we need the static structure function S^(k) of the background and the triplet corre- 
lation function it^ 344 ) with one impurity. Details of the derivation of the above equation are 
given in the Appendix. 

The singularity structure of Eq. ( |2.21| ) is the same as that of the self energy. For real 
frequencies u the imaginary part of (3^j (q) is zero if the energy denominator is negative for 
all values of q. Modes with an energy hu high enough to satisfy the inequality (|2.19|) can 
decay into a phono n-roton mode and the solution has a non-zero imaginary part. 

The long-wavelength limit of the excitation energy defines the hydrodynamic effective 
mass m*r 



2 1.2 



hu 



h z k 



2m* H 



when k — > . 



Inserting this into Eq. (|2.12|) we get 



m 



m 3 



with 



lim 

k^0+ k? 



1 f d*p k.p^ 34 )(p)0p 



(2vr)V 4 ) t(3)(p) + e( 4 )(p) 



(2.23) 



^2.24) 



(2.25) 



where Uq = hk 2 /2m* H . Using Eq. ( |2.23|) we find that in the long- wavelength limit the pole 
strength is inversely proportional to the effective mass, 



lim Z{k) 



m 3 



ni 



(2.26) 



H 



In the so-called "uniform limit approximation"! one neglects all coordinate-space prod- 
ucts of two functions, i.e. we approximate, for example, p^ 34 ^(r , ri)<W 34 )(ro, ri) ~ 
p( 3 )p( 4 )(5u ( - 34 - ) (ro, rj), but convolution products are retained. Then, f3^ {p) has a simple 
form@ 



h^ s^(p) 

P) 2m 3 P S^(p) 



[2.27) 



This together with the equations ( 2.24 ) and ( 2.25|) gives the "un-renormalized effective mass" 
derived by Owenill 
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III. FERMI LIQUID INTERACTIONS 



Up to this point, we have considered only the interaction of one single 3 He atom with the 
host liquid. Further interesting effects arise from the interaction between pairs of 3 He atoms 
and the specific dynamics imposed on the 3 He component by the Pauli principle.0 The most 
obvious manifestations of interactions between the 3 He atoms are magnetic properties^ and 
corrections to the hydrodynamic mass. Both effects provide interesting problems from the 
point of view of low-temperature experiments as well as microscopic many-body theory. 

The dilute mixture is a particularly attractive system for the theorist because many of 
the complicated exchange effects that obscure the theory of pure 3 He are negligible. A new 
effect is that the interaction between 3 He atoms is dominated by the exchange of phonons 
through the host and therefore depends not only on the relative momentum 

between the two particles, but also on the momentum of each individual particle relative to 
the 4 He background. 

Small perturbations of the ground state of a normal, interacting Fermi liquid are de- 
scribed within Landau's Fermi Liquid theoryEi The energy Eq of the system is considered 
to be a functional of the quasiparticle occupation numbers n^ a , and the quasiparticle en- 
ergies are, close to the ground state, given by the variations of that energy with respect to 
n k (J . Thus, the quasiparticle spectrum is 

e (3) (k,.) = ^. (3.1) 

For a symmetric ground state, the e^(k,a) are spin-independent; we shall suppress the 
spin argument whenever possible. 

In the mixture, the single particle spectrum consists of two parts. The first is the 
hydrodynamic interaction of a single impurity atom with the background, leading to the 
"hydrodynamic effective mass" m* H as discussed in the preceding section. The second part is 
due to the interaction with other 3 He atoms. Hence, we can split the quasiparticle spectrum 
(|3.1|) into two pieces, 



e^(k) = e^(k) + e^(k), (3.2) 
Correspondingly, the total effective mass of the 3 He particles is 



h 2 k,F h 2 kp d 



F _ n k f _«_ e (3) , 

m* m* H dk ^ p 



(3.3) 



k=kf 



where kp is the Fermi momentum of the 3 He component. 

Interactions between the 3 He atoms are, to the extent that the 3 He atoms remain close 
to the Fermi surface, described by Landau's Fermi liquid ("quasiparticle") interaction which 
is the second variation 



5e^(k,a) 5 2 E 

ko-,k'cr' ~ 



rvar 
J] 



(3.4) 



n, 



evaluated at the ground state momentum distribution nf^ = 6(k F — k). The quasiparticle 
interaction normally contains a spin- independent and a spin-dependent part, 
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& = £ (AT ± /U ) i (3-6) 



/ko*k'<r' — /k,k' + /k,k' °" * a ' > (3-5) 

which is also frequently written as 

1 

2 

the ex's are Pauli spin matrices. Since the quasiparticle interaction is to be taken for |k| = 
|k'| = hp, it depends only on the angle between k and k', and it is convenient to expand 
that dependence in Legendre polynomials 

& = E ft a) m ■ k') • (3.7) 

i 

The strength of the interaction relative to the kinetic energy is measured by the dimensionless 
quantities 

F £ s(a) = N(0) ff a) = ^f-^ /, s(a) , (3.8) 

where N(0) is the density of states at the Fermi surface. 

At this point, the quasiparticle contribution to the effective mass is normally related to 
the Landau parameter f( or F(. To derive this relationship, one must assume Galilean- 
invariance of the whole liquid. As already observed by Baym and Pethick,0 the argument 
must be treated with caution here because, while the whole mixture is Galilean invariant, 
the properties of the 3 He component are not invariant against motion relative the 4 He back- 
ground. This becomes relevant when one attempts to treat the phonon-mediated processes 
as an effective interaction between the 3 He atoms: As long as this interaction is assumed to 
be static, it can depend only on the properties of the participating 3 He particles. However, 
in reality the interaction is dominated by the dynamic exchange of phonons through the 
background liquid and the motion of the interacting particles relative to that background 
becomes important. The correct definition of the effective mass is therefore the one of Eq. 
( p.3|) ; we will explicitly see the effect of the retardation of the quasiparticle interaction in 
our numerical applications. 

The magnetic response is determined by the change of the quasiparticle energy when a 
magnetic field H is applied. Since the 4 He background does not couple to the spin degrees 
of freedom, the 5e^ 3 ^(k, er) depends, to first order in the magnetic field, only on the spin — 
populations n k (T . Hence the usual relationship of Landau theory, 

&** = Tk fl + ^|ffi/ a )=^| (1 + F a ) , (3.9) 

between the Fq and the spin-susceptibility is applicable. In order to separate the hydrody- 
namic "backflow" effect from the Fermi liquid interaction, we have defined here x* ideal as 
the Pauli susceptibility for 3 He particles with mass m* H . 

A. Variational Fermi Liquid Theory 

Landau's Fermi liquid theory is phenomenological in the sense that it relates Fermi liq- 
uid parameters to physical observables, and derives relationships between different physical 
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observables, but it makes no statements on how to compute these quantities. Methods to 
calculate the quasiparticle excitations and interactions from the underlying Hamiltonian 



( |1.1| ) of the system are provided by microscopic many-body theory. 

At the level of the variational theory, the single particle excitation spectrum is calculated 
by allowing for an occupation of single particle orbitals n^o- in the Slater function $o that 
is different fr 
derivatives hi 
general form 



is different from the Fermion ground state n k . The graphical analysis of the functional 
derivatives has been carried out in Ref. ITS]; the variational single-particle field u(k) has the 



,(3) 



(k,<r) =t®(k) + u(k) + U , (3.10) 



where Uq is a constant determined by the chemical potential, and u(k) a momentum depen- 
dent average field. In the present dilute case, the average field can be written in the form 
of a Hartree-Fock field 

<k) = - —] d 3 rW eS (r) Jo (rk)£(rk F ) (3.11) 

where £(x) = 3 ji{x)/x is the Fermi exchange function familiar from Hartree-Fock theory, 
and Weff(r) is a local effective interaction which can be constructed diagrammatically by the 
summation of (F)HNC diagrams in the mixture.i'0 For further reference, we also formulate 
the average field in momentum space, 

/d 3 o 
(27r)V 3 ) ra lq- k l- g Wgffta) » ( 3 - 12 ) 

where W e s(q) is the dimensionless Fourier transform W e g(q) = J d 3 rW e s(r) exp(z k ■ r). 

In the same approximation, the quasiparticle interaction can also be expressed, up to an 
additive constant, in terms of the effective interaction W e g{r): 

/k™ k v = ^ [const. - W eS (\k - k'|) 

= fLy + fl V * ■ <*> . (3.13) 

The constant first term in Eq. ( |3.13| ) comes from the variation of Uq with respect to the 
occupation number n^y, it is related to the compressibility of the 3 He component and is 
therefore not our concern now. 



B. Correlated Basis Functions 

It has been known for quite some timeE! that the naive implementation of the FHNC/EL 
theory (or, for that matter, any variational wave function of the type ( p..3| )) to the problem 
of pure 3 He has a number of severe problems. One is that it leads to an effective mass 
of pure 3 He that is less than one@ which is in sharp contrast to experiments.0@ Another 
problem is that it predicts an instability against spontaneous spin-polarization. The physics 
of mixtures is, of course, entirely different from that of a pure 3 He liquid. But we will see 
that the technical approximations implicit to the Jastrow-Feenberg wave function ( [L3|) lead 
to practically the same problems in the mixture. 
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The deficiencies of the variational theory are, as we shall see more explicitly below, due 
to the fact that the wave function (|1.3Q describes the average correlations between particles 
quite well, whereas it is insensitive to the specifics of the correlations in the vicinity of the 
Fermi surface. The cure for the problem is known to be Correlated-Basis Functions (CBF) 
theory.i A complete, non-orthogonal basis of the Hilbert space is built on the ground state 
wave function (|1.3|) by defining 

M{rS a) }) = e^ Ca) »$ m ({rf >}) , (3.14) 

where the |$ m ({rj 3 ' ) })| is a complete set of Slater determinants for the 3 He component. 

Within the basis |^m({ r i })}> perturbative expansions can be derived in much the same 
way as in conventional Rayleigh-Schrodinger perturbation theory, and diagrammatic struc- 
tures like ring-, ladder-, or self-energy diagrams can be defined and summed to infinite order. 
Note that perturbative corrections are applied to the Fermion components only; this is le- 
gitimate since the Jastrow-Feenberg wave function is in principle exact for Bosons, whereas 
it is not exact for Fermions. 

The general strategy of CBF theory has been described in detail in many places;! the 
basic results needed for our analysis have been derived in Ref. [23| Without going into 
the details of that theory, we assert that there is a one-to-one correspondence between 
the compound-diagrammatic elements of the FHNC theory and that of ordinary, Green's 
functions based perturbation theory. Loosely speaking, CBF theory can be interpreted as 
a microscopic procedure for the calculation of static effective interactions. For our purpose, 
the application of CBF theory can be summarized simply by stating that, in the evaluation 
of the energy, the chain- diagrams of FHNC theory are replaced by the ring-diagrams of the 
RPA theory, using a diagrammatically defined static effective interaction which we identify 
with the "particle-hole interaction". 

Since we are predominantly interested in variations of the energy, we can formulate the 
theory in terms of a Green's functions approach where the effective interactions come from 
the FHNC analysis, examine the relationships between the FHNC and the CBF approxima- 
tion, and compare the results of different theories. 

In perturbative theories, single particle properties are described by a complex self-energy 
H(k,u); the single particle spectrum e(k) is obtained from the solution of the equation 

e(k) = t (3 \k) + Z{k,e{k)/h) . (3.15) 

If only single phonon coupling processes are considered, the self-energy H(k,E) is given in 
the so-called GOW-approximationBi! 

m U )=iJ ^0^- G (0) (|k - q| , U - u/) V cS (q, u>) (3.16) 



where 



G^(k,u) = - 1 , + !^ CT (3 1?) 

v ' hu -t®(k) +ir) hu -t( 3 )(fc) -in v ' 



is the free single-particle Green's function and 
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V eS (q,uj) 



a/3 



(3.1* 



is the effective, energy dependent 3 He- 3 He interaction. In Eq. ( |3.18| ), the Vp_ b (q) are the 
local, particle-hole irreducible interactions, and {.Qi u ) ls the density-density response 
matrix. In fact, the GOW approximation is simply the variation of an RPA energy expec- 
tation value with respect to a single particle Green's function. The key quantities V^f h (q) 
can in principle be defined diagrammatically in Green's functions based theories, but it is 
impractical to calculate these interactions without further approximations. CBF perturba- 
tion theory solves this problem by giving unambiguous prescriptions^ for the calculation of 
static approximations for these quantities. 

To separate the "hydrodynamic" and the "Fermionic" component of the self-energy, we 
rewrite the single-particle Green's function as 



G<°>(Jfe, 



1 



%u-t$)(k) +ir] 



1 



+ 



1 



%u - t^ (k) +irj %lo- t( 3 ) (k) - irj 



and the self-energy in the form 



d 3 qd(hu>') ^(o) 



E H (k,uJ 



(27r) 4 p (3) 
d 3 qd{hu;') 
(2vr)V 3 ) 



G^(\k-q\,u-uj')V cS (q,uj'), 
GP(\k-q\,u-u')V cS (q,Lj') 



(3.19) 

(3.20) 
(3.21) 

(3.22) 



The "hydrodynamic" part ( |3.21| ) of the self-energy leads, in the dilute limit — > 0, to an 
effective mass that is identical^ to the one obtained in the previous section in the "uniform 
limit approximation" derived in Eqs. ( [2.24 ) and ( |2.27| ). This approximation is quantitatively 
not quite sufficient, but it contains, at a semi-quantitative level, the basic physics of impurity 
motion, in particular the coupling of the impurity to the excitations of the host liquid. The 
method of time-dependent correlation functions described in Sec. [J introduces a systematic 
procedure to improve upon this approximation. 

Let us now focus on the Fermion term, Eq. ( |3.22| ). The energy integration yields the 
compact form 



E F (k,u) = - 



d 3 q 



n\q-kiaV cS (q,u - t (3) (k - q)/h) 



(3.23) 



(2vr) 3 p( 3 ) 

To calculate the spin-susceptibility we can directly apply Eqs. ( |3.8| ) and ( ^.9[ ) and arrive 



at 



n/5 



2p( 3 ) 



d9 



sin(0)y eff (|k F -k' F |,O). 



(3.24) 



In principle, we would also need to calculate the variation of the kinetic energy t^(k — q) 
with respect to the quasiparticle occupation number n^o-, but in the dilute limit of the mix- 
ture, this quantity is dominated by the hydrodynamic backflow relative to the background 



11 



and can, to good approximation, be replaced by a free single-particle spectrum with an 
effective mass m* H . 

The comparison of Eqs. ( |3.24| ) and ( |3.23j ) demonstrates the difference between the Lan- 
dau effective mass that would be calculated from the I = 1 equivalent to Eq. ( |3.24|) , and the 
momentum derivative de(k, a) / dk\ k=kp . In the latter case, one also includes the momentum 
dependence of the energy appearing in the effective interaction V e g(q,tu — t^(k — q)/h). 

We have now reached the point where the quasiparticle interaction is expressed in terms 
of a local, but energy dependent effective interaction very similar to the FHNC theory, cf. 
Eq. ( |3.13| ). The relationship to Eq. ( |3.12| ) is now apparent: Hpik^u) could be rewritten 
as a Hartree-Fock expression of the form (|3.12|) if the effective interaction V e ff(q,u>) were 
energy independent. The question that arises naturally is what the relationship between 
the energy-independent effective interaction W c g(q) and the energy dependent interaction 
V e g(q, uj) is. We will address this question in the next subsection. 



C. Connection between HNC and CBF Theories 



We here outline the general manipulations that lead from the compound-diagrammatic 
elements of the (F)HNC theory to the Green's functions expressions. The construction 
of energy-independent, local effective interactions is one of the key steps in the Feynman- 
diagram based parquet-diagram theory@i3 that is needed to show the equivalence to the 
HNC/EL theory, and we shall demonstrate here that the relationship between (F)HNC and 
parquet theory reach actually farther than expected. 

For simplicity, we restrict our derivation to the case of two interacting 3 He impurities 
within the 4 He host liquid. In this case, we need to retain only the 44-channel of the response 
function. In this limit, the effective interaction is 



V cS (q,u) = V^(q) + V£S( 9 ) X m (q^) V<$(q) 



vS(q)+V^>(q) 



(34), 



2t«(g) 



(M 2 - (e (4) (<?)) 5 



(3.25) 



The prescription from parquet-theory to make this energy-dependent interaction local is as 
follows: Construct the RPA static structure function 



4p1(<z) 



d(hu) 



2tt 



3m 



Xo 33) (?, <*>) + Xo M> (q, w) V eS (q, w)Xo^ j (9, w 



,(33) 



(33), 



1 - 



where 



f0%) 



xo 3 \q^) 



V$$(q)-2 



(t®(q) + eW(q)y 



(3.26) 



2t^{q) 



(hcu) 2 -(t®(q)y 



(3.27) 



is the response function of the non-interacting 3 He gas. Also, construct the ladder approx- 
imation for the same quantity in terms of a different and yet unspecified local effective 
interaction, say Vl(<?) 
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o(33) 



ladder 



(?) 



d(H(j) 
2vr 



3m 



#)(g) 



(3.28) 



Now choose an average frequency u)(q) such that these two forms of the static structure 
function are identical for 



V L (q) = V eS (q,u) 



One obtains 



(hu>) 



eW(g) (t^nq) 
eW(g) + 2t( 3 )(g) 



(3.29) 



^3.30) 



This is the generalization of the expression given in Ref. |2?| to a dilute 3 He- 4 He mixture. 
Inserting fiu into the energy- dependent effective interaction leads directly to the FHNC ap- 
proximation ( |3.11| ), in other words 



W eS (q) = V e{i (q,ij) = V L (q) 



(3.31) 



This final part is of a somewhat technical nature and requires the full (F)HNC formalism, 
in particular its use to calculate effective interactions for use in CBF theory!! and the 
connections with the optimization problem; it is therefore skipped here. 

The equivalence (|3.31 ) is the second essential theoretical result of this paper. The result 
is remarkable in the following sense: The concept of an "average energy" has been introduced 



in Refs. |26] and |27| for a one-component Bose liquid, and for finding a local approximation for 
the chain diagrams for the purpose of summing ladder diagrams. We find here - as in Ref. |2| 
- that the very same concept applies also in a different situation, namely in demonstrating 
the relationship between the quasiparticle interaction in FHNC theory and in the RPA. 
We stress that this result is an observation on how the static FHNC approximation and 
the GOW approximation of the self-energy are related. It does not mean that the FHNC 
approximation, which has been derived from looking at the static structure function, and 
which is capable of reproducing the static structure function with very good accuracy, is 
also adequate for the calculation of single-particle properties and Fermi liquid effects. 



IV. APPLICATIONS 

A. Hydrodynamic Effective Mass 

The first task is the calculation of the hydrodynamic effective mass. This quantity de- 
scribes the interaction of a single impurity atom with the 4 He background. Due to the high 
density of the background, the GOW approximation recovers only about 80 percent of the 
effect and the more complete treatment of the time-dependent pair correlations described 
in Sec. [11] improves this result considerably. 

One of the formal results of our analysis is that the equations (|2.21|) , ( |2.13| ) and (|2.22|) 
treat the phonon-roton spectrum as a Feynman-type spectrum, e^(k) = h 2 k 2 /[2rri4 i S(k)], 
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which is known to lie by a factor of two too high in the roton region. To include the 
backflow correction also into the motion of the 4 He particles one should add fluctuating 
triplet correlations into the wave function ( |2.3| ) and solve the triplet continuity equation. 
This leads to a rather complicated integral equation with uncertainty in the treatment of the 
four-particle distribution function. The hierarchy of continuity equations must be truncated 
at some tractable level. A reasonable shortcut is to modify the integral equation (|2.21| ) and 
the self energy ( |2.13| ) by using the experimental phonon-roton spectrum in e^(k), replacing 



the impurity mass by its effective mass in the impurity kinetic energy term, t^ 3 \k), and then 
ignoring the higher order correlations. For the triplet distribution function we have used 
the convolution approximation together with the triplet correlation function as described in 
the Appendix!! The self-consistent solution of Eq. (|2.24j ) leads to a rather good prediction 
of the hydrodynamic effective mass over the whole density regime that is experimentally 
accessible. The results are shown in Fig. || together with the fits to the measurements by 
Yorozu et a/.0 and Simons and Muellei@. More details of those results are shown in Table 
|, which will be discussed in the next section. 

Next we need to determine the momentum dependence of the hydrodynamic effective 
mass; it is necessary to study this effect because at finite concentration the particles at the 
Fermi surface have a finite momentum. For that purpose, we have calculated the disper- 
sion relation from Eq. ( |2.12| ) and have determined the "momentum dependence" of the 
hydrodynamic mass by writing the spectrum in the form 



In Fig. ^| we plot the quantity 



1 



k 



m* H (0) 



- 1 



(4.2) 



which turns out to be constant for a wide momentum range implying quadratic momentum 
dependence of the effective mass. Our value of this constant is 0.114 at zero pressure. It 
has also been estimated experimentally by Fak et aZ.S They find the value 0.114±0.01 at 
the pressure P=0.1 atm and 1% concentration by fitting the momentum dependence of the 
spectrum in the range 0.9A- 1 < k < 1.65 A" 1 . This is in very good agreement with our 
result. The density dependence of this constant is weak both in our calculations and in 
experimental analysis. A slightly larger value 0.19 was obtained recently by Fabrocini and 
PoiiM That is due to the fact that their calculated spectrum is slightly softer and ours. 

Different_parametrizations of the experimental spectra for k < 1.7 A have been 
suggested.!!'!! Fig. ^ shows a comparison between our dispersion relation and these 
parametrizations. Our spectrum is only slightly stiffer, but the agreement appears to be 
satisfactory for all practical purposes. 



B. Damping of Impurity Motion 

At higher momenta and associated energy, it becomes kinematically possible that the 
3 He decays into a roton and a lower energy impurity mode. At this point, the self energy 
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becomes complex and no solution for the pole equation Q2.12j ) can be found.il The dynamic 
structure function has no longer a delta function peak; instead the energy is distributed 
over a range of modes. Therefore for momenta k > 1.7 A we plot in Fig. |3| the maximum of 
the dynamic structure function, which shows that the impurity mode crosses the phonon- 
roton spectrum smoothly. The peak broadens substantially before the second crossing at 
k » 2.25A. 

The pole strength defined in Eq.fl2.20p is calculated by numerically differentiating the 



self energy. That is very reliable because the self energy is a smooth function of & in 
that frequency range. The results are plotted in Fig. [|. As pointed out in Eq.( |2.2(j| ) its 
long wavelength limit is inversely proportional to the effective mass. The monotonically 
decreasing behavior with decreasing wavelength is partly due to the increase of the effective 
mass, also shown in the figure, and partly due to the increasing contribution of the phonon- 
impurity channels at higher frequences. Fak et al. have measured the area under the 
particle-hole peak of the 1% concentration mixture for the momentum range 0.9A -1 < 
k <1.65A _1 . These results shown in Fig.|] are in very good agreement with our calculations. 
The measurements with 5% concentration suggest only a weak concentration dependence 
and thus the comparison with the zero concentration limit calculations is justified. 

In Fig. [5] we show the grey scale plot of the low energy part of the S(k, us) for u < 25 K. 
The pole contribution of the elementary excitation mode has been artificially broadened 
by convolution with Gaussians in such a way that the area under the peak is equal to the 
strength of the pole. In the figure we also show the experimental phonon-roton spectrum 
and the continuum boundary. The elementary excitation mode crosses that boundary at 
1.72 A -1 and the decay into a roton becomes possible. The dominant peak in S(k } u) crosses 
the phonon-roton spectrum and looses its strength to higher lying modes. The mode slightly 
above 20 K is clearly visible in the figure. 

Damping of the impurity motion in the roton region is more clearly shown in Fig. ||. 
There we plot S(k, u) as a function of uj for some typical values of k = 1.9, 2.1 and 2.3 A -1 . 
The peak, which is a delta function when k <1.72 A -1 broadens and looses its strength when 
the roton minimum is past. In Table || we give the half widths of those peaks. From them 
we can estimate the life time A of the mode. The mode propagates with the group velocity 
and from that we can calculate the distance A the impurity can travel before it looses its 
energy to a roton, 

h 2 k . . 

A = 77T^ • 4.3 

The propagation distance is typically a few Angstroms as also shown in Table |T|. 

The penetration depth of a 3 He atom at energies that are high enough such that the 
impurity atom can couple to the roton are important for the interpretation ofrecent mea- 
surements of the momentum transfer of single 3 He atoms from 4 He clusters^. In these 
experiments, a significant increase of the momentum transfer has been observed that is 
consistent with our energies. The very rapid drop of the penetration depth shows that prac- 
tically all 3 He atoms colliding with a 4 He cluster will be absorbed from clusters as small 
as 100 atoms, which have a diameteiS of about 25 A. One of the reasons for this is the 
comparison of the cluster size with the penetration depth, the other is that the excitation 
spectrum of even such small clusters is remarkably similar to that of the bulk liquid. The 
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argument does not yet include the inelastic coupling of the impurity atom to the lower-lying 
surface excitations, which can be siza ble§. 



C. Fermi Liquid Effective Mass 

We now turn to our microscopic calculation of Fermi liquid effects. The necessary ingre- 
dients of the theory — the effective interactions V^°_^{q) and W c g(q) as well as the Feynman 
spectrum e^(k) — have been obtained in Ref. |2|. There are two sets of accurate experi- 
mental data, those of Yorozu et a/.il and those of Simons and Mueller.il These experiments 
differ by the pressure normalization, but lead otherwise to very similar results. 

Clearly, the effective mass is dominated by the hydrodynamic backflow as calculated 
above. To eliminate any uncertainty caused by inaccuracies in that calculation, we have, for 
the further calculations, not used the theoretical values obtained in Section || and shown in 
Figs. and ^, but rather considered the zero-concentration limit of the hydrodynamic mass 
as a phenomenological input. After the concentration dependence was calculated from the 
Fermi liquid contributions as described below, we have fitted our theoretical values to the 
experiments of Refs. |28| and ^ to optimize the overall agreement and then calculated the 
zero-concentration extrapolation. 

Three calculations have been carried out to determine the Fermi liquid contributions to 
the effective mass of the 3 He component as a function of concentration and pressure. The first 
calculation applied the simple FHNC/EL theory and the static effective interaction ( p,12|) . 



To account for the hydrodynamic backflow, one must supplement the Fermion contribution 
( p,10| ) by the hydrodynamic contribution hunik)] then the spectrum has the form 



e( 3) (k) = hu H {k) + u{k) + U , (4.4) 



where the Fermi correction u(k) is given in Eqs. ( |3.11| ), ( |3 . 1 2| ) . When treated this way, the 



concentration dependence of the effective mass derived from the spectrum (f4.4|) is visibly 
steeper than the experimental one, as seen in Figs. [7] and [8]. 

In the next step, we calculate the effective mass using Ks(k, 0) as quasiparticle- 
interaction equivalently, setting Q = in Eq. (|3.31|). This form of the self-energy relaxes 



the approximations made by the FHNC theory since it takes the effective interaction at the 
Fermi surface and not at an average energy. We see in Figs. |7| and [5] that the agreement with 
the experiment is indeed improved; the approximation recovers about half of the discrepancy 
between the FHNC approximation and experiments. 

Finally, we have carried out self-consistent calculations of the effective mass. It is suffi- 
cient for that purpose to assume a single-particle spectrum of the form t^ 3 \k) = h 2 k 2 /2m* 
in the Green's function ( |3.17| ) and, consequently, in Eq. ( |3.23| ); note that the hydrody- 
namic mass is included in the Green's function. This effective mass is then determined 
self-consistently by requiring that the spectrum e^(k) determined by 

e^{k) = hu H {k) + ll F [k^)i (4.5) 

can be fitted by the same effective mass that has been used in the self-energy. This calcu- 
lation provides a very good agreement with the experimental data as seen in Figs. [7| and || 
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The agreement is worst for the pressure 10 atm and the data of Ref. but we note that 
there is a non-monotonic behavior of the slope of the data as a function of pressure, and it 
might be interesting to re-examine this pressure regime experimentally. 

Carrying out this self-consistent procedure, we arrived at the following interpolation 
formulas for the hydrodynamic mass: 

^) = 2.18 + 2.43 r + 2.67 r 2 - 1.17 r 3 (4.6) 

1^3 J expt 



for the data of Ref. and 

m". 



H \ rt 1 r i o i c„ i a A<-! Ji 



2.15 + 2. 16r +4.47r 2 (4.7) 

m 3 / expt 

from those of Ref. [29. Here, r = p^ /pq — 1, p^' is the 4 He density and po=0.02183A~ 3 is its 



value at the saturation vapor pressure. Typically, the discrepancy between the two different 
extrapolations is 0.03. These extrapolated hydrodynamic masses are shown, together with 
our theoretical calculation of Sec. [II], in Fig. [|. The theoretical values are throughout the 
full density regime about 0.05 below the experiments which is quite satisfactory given the 
level of approximations. 

To produce Figs. |7| and |] we have used - as stated before - the hydrodynamic mass 
given in Eqs. fl4.6|) and ([4.7|) . Since the calculations were done for fixed densities at each 



concentration and the experiments were done for fixed pressure, we have used the experi- 
mental pressure-density relation of Ref. [l| to make the conversion. Our calculations predict, 
at low concentration, a visible curvature of the effectivemass as a function of concentration, 
similar to the predictions of Bashkin and Meyerovich.tJ Hence, we are not convinced that 
linear extrapolations are a legitimate means to determine the hydrodynamic mass unless 
concentrations are significantly lower than those examined in Ref. Such a curvature is 
implicit to the Fermi functions. Already the simple approximation ( f4.4|) would lead to a 
behavior 



m (XI 



m* H + ax 2/3 + bx + cx 5/2 + dx 7/3 .... (4.8) 



The numerical coefficients a . . . d can be calculated from the moments of the potential, but 
such an expansion provides valid results only for very small concentrations and a global fit 
of the form fl4.8|) to the calculated data is more accurate. In Table |HI| we list their values 
for different pressures obtained from the least square fit to the fully self-consistent solution 
ofEq. Q. 

On the other hand, we see no drop in the effective mass at stable concentrations as, 
for example, proposed by Hsu and Pinesll. The monotonic behavior is conclusively shown 
experimentally by Simons and Mueller^ The slight downward curvature of the effective 
mass indicates that such a drop will happen eventually; again the same conclusion can be 
drawn from the analytic form (|4.4j ), but it does not happen at the experimentally accessible 
stable concentrations. 
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D. Magnetic Susceptibility 



We base our analysis on the susceptibility measurements by Ahonen et a/..E3 The spin 
susceptibility of the 3 He component is determined by both the effective mass and the quasi- 
particle interaction in the spin channel,0@ cf. Eq. (|3.9| ), see also Ref. |14] for further 
discussion. To extract the antisymmetric Landau parameter Fq from the spin susceptibility, 
the effective mass must therefore be known from an independent measurement. Again, one 
finds that the spin susceptibility is vastly dominated by the one of the free Fermi gas of 
particles with an effective mass m* H , and one needs very accurate measurements to extract 
information on the quasiparticle interaction. Ref. has used the best values for the effective 
mass available at that timeEil to calculate Fq from Eq. ( |3.9|) . 

To turn to the theoretical description, let us first look at the full magnetic susceptibility 
obtained from Eq. ( |3.9| ) and the effective interactions determined from our microscopic 
theory. Parallel to the effective mass calculation, we use both the static effective potential 
W e s(q) and the RPA effective interaction V e g(q, 0) to calculate the spin susceptibility. The 
experimental data were given on an arbitrary scale, we have therefore scaled each of our sets 
of results with a global factor to provide the best overall fit to the data for the concentrations 
0.27% and 1.33%; the results are shown in Fig. || For comparison, we also show the 
susceptibilities that one obtains from a free Fermi gas approximation with the effective 
mass calculated in the preceding section. Up to five percent concentration, the results look 
satisfactory; the theoretical results from the dynamic theory and those from the free Fermi 
gas approximation basically bracket the experimental data. There are some deviations at 
the phase separation concentration; here the free Fermi gas model is somewhat worse than 
the dynamic interaction. These deviations can, on the theoretical side, be caused by the 
fact that the convergence of our FHNC approximations becomes worse at higher Fermion 
concentration. Part of the disagreement is also because the data at lower pressures are 
given the phase-separation concentration and only those at higher pressures at 8.8 %. The 
limiting solubility at zero pressure is approximately x = and increases until it reaches 

a maximum of about x = 0.095 near 10 atm. The theoretical calculations were, on the other 
hand, all done at 8.8 % concentration. 

A remarkable result is that the susceptibilities derived from the static interaction W e s(q) 
are clearly poor. In fact, the independent particle approximation provides better agreement 
with experiments than the static variational theory. This is not entirely unexpected in 
view of the above analysis of the effective mass, and also in view of the much more severe 
deficiencies of the Jastrow-Feenberg wave function for magnetic properties of pure 3 He.0 

The purpose of the points made above was to raise a general awareness of the difficulty 
of obtaining information beyond the effective mass from susceptibility data in mixtures; this 
difficulty has already been pointed out by the Ahonen et ai; further concerns on the preci- 
sion of the pressure normalization of that work were raised by Rodriges et a/.0. To extract 
the Landau parameter Fq from these data, we have followed the procedure of Ahonen et al, 
but also used the more recent effective masses data of Refs. and ^ discussed above. The 
resulting Landau parameter Fq is shown, together with the original data of Ref. |39], those 
obtained from the raw data with the effective mass fit (|4.7|) , and our theoretical results, in 
Fig. [10]. Obviously, the conclusions one can draw from the same measurements depend visi- 
bly on what effective mass is used to extract Fq from the spin-susceptibility. The agreement 
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between theory and experiment is significantly improved by using the more recent effective 
mass data. There is still some vertical offset but the concentration dependence of our theo- 
retical results is in fact quite good; the theory predicts a somewhat stronger concentration 
dependence of Fq than is seen experimentally. Most of the remaining disagreement can be 
removed by using a slightly different hydrodynamic effective mass. To demonstrate this, we 
have repeated the procedure of the preceding section and fitted the hydrodynamic effective 
mass such that, at finite concentrations, the theoretical Landau parameter is reproduced. 
The best fit to the data is reached by choosing the hydrodynamic effective mass ratios 
given in the last column of Table |. It appears that, apart from a slightly larger curvature 
suggested by theory, there is little room for improvement. It is clear that any attempt to 
extrapolate an effective mass from available susceptibility data is uncertain because there 
are simply not enough data available to carry out such extrapolations with confidence. 

In this connection we also would like to stress the rather strong concentration dependence 
in the low-concentration regime. It indicates that, for the purpose of extracting Fq from 
susceptibility measurements, even a concentration of 0.27% is far from the dilute limit. Note 
that — if one assumes no readjustment of the hydrodynamic mass — the concentration 
dependence of Fq as x — > must be assumed to be even stronger. Similar to the effective 
mass, the concentration dependence of the Landau parameter Fq is most clearly discussed 
by writing Eq. (|3.24|) in coordinate space: 

F o=~^2 fd 3 rV eS (r,0)f (rk F ). (4.9) 



27r 2 h 2 

One source of the concentration dependence of F$ is kinematic dependence coming from 
the explicit appearance of the Fermi wave number kp in Eq. fl4.9|) , which is caused by the 
Pauli principle, as well the effective mass ratio in front of the integral. The second is the 
implicit dependence of the effective interaction V e g(r, 0) on the concentration. We found, 
however, that this dependence is negligible within the experimentally accessible regime under 
consideration here, and that V c s(r, 0) is quite well represented by the interaction between a 
single pair 3 He atoms within the host liquid. In other words, the concentration dependence 
is almost entirely due to the kinematics dictated by the Pauli principle. 

The coordinate-space representation ( [4.9| ) suggests a concentration expansion similar to 
Eq. (|4.8|) for the antisymmetric Landau parameter. Since Eq. (|4.9| ) suggests a natural 
factorization into an effective mass ratio and an interaction term, we expand 



F?(x) = a x 1 '* + bx + c x 5 / 3 + d x 7 / 3 . (4.10) 

m 

The density-dependent parameters entering this fit are given in Table [TV]. From the 



form of our results shown in Fig. |10| it appears that the first two coefficients should suffice; 
inclusion of two more terms improves the fit slightly. In principle, one can obtain the first 
two coefficients from moments of the effective interaction. In practice, however, only the 
x i/3 t erm gi ves a faithful approximation of our results up to 1.33 % concentration. 

Effective quasiparticle-interactions have in the past been discussed and parameterized 
mostly in momentum space. This seems to be appropriate because the quasiparticle in- 
teraction is indeed a low-momentum property of the effective interaction. On the other hand, 
the properties of these interactions are more intuitively described in coordinate space. The 
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effects that contribute to the effective interaction have been discussed by Aldrich, Hsu, and 
Pines;@@ these are (a) core exclusion due to short-range correlations, (b) an increase of the 
core-size due to the kinetic energy needed to bend the relative wave function to zero within 
the core, and (c) an increased attraction due to the presence of other particles in the area 
where the potential is most attractive. We show the effective interactions entering Eq. Q4.9J1 
as a function of density at 5 % concentration in Fig. [IT]; the underlying bare Aziz-potentiaB 
is also shown for reference. All the effects predicted by Aldrich and Pines are clearly seen 
in Fig. [n] We also note that these overall features of the potential are quite resilient; the 
density dependence of the interaction is relatively weak and the concentration dependence 
practically invisible on the scale of Fig |TT[ On the other hand, the differences between the 
"dynamic" and the "static" (FHNC) approximations are larger than the variation of the 
potentials over the experimentally accessible density regime of interest here. 

Finally, a word is in order concerning the behavior of the dynamic interaction within the 
core region. V e R (r, 0) should be zero in this regime. In practice, it does not vanish vanish, this 
is due to the "RPA"-like approximation ( 3.25 ) for the induced interaction. To be completely 
consistent, one would have to solve the ring- and ladder diagrams self-consistently using 
the full Fermion propagators, and not the "collective" approximation. This task has not 
been accomplished yet. We have therefore set the interaction V e g(r,0) to zero within the 
core region. Doing so or not has no visible consequence, which lends credibility to our CBF 
treatment that introduces the correct particle-hole propagators a posteriori in a perturbative 
way. 



E. Scattering Matrix and Phase-Shifts 

To study the transport phenomena and potential superfluidity of the 3 He atoms in the 
medium, one needs to consider the scattering of two 3 He atoms at momenta k and k'. 
Again, we must consider both the hydrodynamic backflow, and the direct interaction which 
is dominated, for long distances, by phonon exchange and for short distances by the bare 
interaction between particles.0@ A series of experiments to find a superfluid phase transition 
of the 3 He componentEHEll exhibit no evidence of superfluidity down to a temperature of 
about 100 /iK; theoretical estimates@@ of the transition temperature range over several 
orders of magnitude. 

Our calculations are similar to those of Owenjll but use the more accurate ground-state 
results of Ref. The interaction entering the scattering equation is not the same as the 
one used above for the calculation of Fermi liquid parameters. The reason for this is that 
in the above calculation, short-range correlations are included by calculating the "parallel- 
connected" diagrams or, in the language of perturbation theory, the "ladder diagrams". In 
the scattering equation, these short-range correlations are dealt with by solving an effective 
Schrodinger equation which reduces, for the case that both scattering particles are in the 
many-bo dvground state, to the Euler equation of the ground state theory. The scattering 



equation isc 



where 



h 2 

V 2 + V; cat (r)-£ 0(r) = O, (4.11) 

m 3 
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V scat (r) = V(r) + AV e (r) + wi(r) . (4.12) 

Here, V(r) is the bare potential, wi(r) the "induced interaction" of the FHNC theory, and 

AV c (r) a correction due to elementary diagrams and triplet correlations. In the dilute limit, 
the pair-distribution function g 33 (r) between pairs of 3 He atoms is given by the zero-energy 
equation 

h 2 

V 2 + V scat (r) 



m 3 



<? 33 (r) = 0. (4.13) 



Eq. (|4.13|) is one of the central equations of the HNC/FHNC theory, see, for example. Ref. 
for derivations and discussion. 

The only difference to Owen's work is so far that we have included elementary dia- 
grams and triplet correlations in the ground state, and hence the induced potential wi(r) 
also changes. The effective potential is shown, for three representative densities, in Fig. 
|T2"| . Again, this effective interaction is almost independent of the 3 He concentration; the 
quasiparticle scattering potential depends therefore on the density and not on the concen- 
tration. Although technically not necessary, it would be acceptable to use (as Owen) the 
low-concentration limit of the effective interaction to calculate scattering phase shifts. 

One might now be led to argue that the wi(r), which describes phonon exchange, should 
also be a dynamic, energy-dependent interaction as discussed above. This is in principle cor- 
rect, but unfortunately there is presently no practical way include such dynamic effects. The 



reason for this is that the ground state equation Q4.13j ) states that J g 33 (r) is a zero energy so 



lution of the scattering equation ( 4.11 ). Changing the induced interaction a posteriori leads 



to an inconsistent low-energy behavior of the solution and to spurious bound states. Hence, 
a better calculation must await the developments of a complete parquet-diagram theory that 
includes energy-dependent interaction all the way through the diagram summations. 

Returning to the original problem, we have solved the scattering equation ( [4.1 1| ) as a 



function of energy in the I = 0,1, and £ = 2 channels. The corresponding phase shifts Sg are 
shown, as a function of density, in Fig. |H| 

The definition of the scattering amplitude isfLZl 

T t {k) e-*< = -^-(21 + 1) sin(fc) , (4.14) 

777-3 /C 

where k = y/2m 3 E/h is the relative momentum. 

The most interesting application of our results is the estimate of a potential superfluid 
phase transition. Manifestly microscopic many-body theory is still at a rather unsatisfactory 
state when it comes to predicting a superfluid phase transition in 3 He, and the mixture 
problem is no exception, possibly because of the lack of a self-consistent parquet-like theory. 
At this point, we can only use the scattering phase shifts in a weak-coupling approximation 
to estimate the critical temperature of the phase transition. 



k B T c m Ep exp 



2*2 



n 2 h 



m* k F T e (2E F ) 



(4.15) 



h 2 k 2 o 

where Ep = ^ — f is the Fermi energy and Tg is evaluated at 2Ep because two He with Fermi 



energy Ep are forming a Cooper pair. Eq. (|4.15|) gives only a rough estimate of the critical 
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temperature since it overestimates the critical temperature by two orders of magnitude in 
pure 3 He. 

Our results for the critical temperature are shown in Fig. ^ for two representative 
densities. The general trend is that the critical temperature for p-wave pairing increases 
strongly with concentration, while the temperature for the s-wave decreases; this is mainly 
due to the explicit appearance of the Fermi momentum in Eq. ( |4.15|) and the resulting 
dependence of the scattering matrix elements Tg(E). T c decreases with increasing pressure 
which can also be seen from the decrease of T shown Fig. [L4| as function of the energy. 

Our results are an order of magnitude smaller than Owen's which we have reproduced, 
but significantly larger than the values proposed in Ref. |39|. The difference to Owen's results 
is mainly due to the smaller scattering matrix elements in both s and p wave channels due 
to the more quantitative treatment of many-body correlations. However, the wide spread 
of results demonstrates drastically the very sensitive dependence of the critical temperature 
on the matrix elements. Therefore, while we believe that the scenario shown in Fig. 15 is 
qualitatively accurate, we would like to be cautious about the quantitative validity. 

This is consistent with our calculation; however we prefer to be cautious to call the 
results of our rather simple calculation quantitative for reasons explained above. 



V. SUMMARY 

We have in this paper analyzed various procedures for calculating the quasiparticle inter- 
action of the 3 He component in a 3 He- 4 He mixture. The technical parts of the calculation 
were based on the equation of motions method, optimized (F)HNC theory of mixtures and 
on CBF theory to infinite order. The close relationship between the theories has been 
discussed. 

Our results for the single-impurity spectrum are quite satisfactory. One might object 
at this point that our use of the experimental phonon-roton spectrum leads us away from 
manifestly microscopic many-body theory, but we feel that such shortcuts are legitimate to 
eliminate tedious and unrewarding computations whose outcome is basically known. 

Turning to the Fermion aspect of our calculations, we have seen that the naive FHNC 
theory is generally unsatisfactory for the prediction of Fermi liquid effect compared with 
level of accuracy that has been obtained for the microscopic calculation of ground state 
properties of simple quantum liquids and quantum liquid mixtures. We have also noted a 
systematic, stepwise improvement of the theory when dynamic and Fermi surface specific 
effects are included. 

The situation is clearest for the magnetic susceptibility where the FHNC approximation 
takes the effective interaction at the energy hu(k) whereas it should have been taken at 
zero energy. We stress that this is not a shortcoming of the (F)HNC theory which is a 
specific method to sum large classes of diagrams. It is a problem of the Jastrow-Feenberg 
function per se. It is also clear how infinite order CBF theory resolves the problem in a 
natural and elegant way, whereas it would be quite difficult to describe the same physics by 
attempting to construct "better" variational wave functions for the ground state. The same 
effect is also the explanation for both the miserable^ performance of simple variational wave 
functions for the effective mass in pure 3 He, and for why magnetic properties of pure 3 He 
are not well described by the naive application of variational wave functions. The power 
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of the CBF approach lies in the possibility of combining two generic many-body methods: 
Green's function concepts are used for the examination of subtle, energy dependent effects 
and variational methods for the unambiguous determination of static effective interactions 
whenever such interactions are appropriate. 

In conclusion, it appears that one understand the concentration dependence of the two 
Fermi liquid parameters Fq and reasonably well with a local quasiparticle interaction 
if retardation effects are included. The remaining discrepancy between our results for the 
effective mass and the magnetic susceptibility seems to prevail at even very low concentra- 
tions where Fermi liquid effects are particularly well described by our theory. The situation 
becomes even more disturbing considering the fact that the data at 0.27% and 1.33% con- 
centration seem to indicate a smaller curvature than theory, which should then turn steeper 
in order to reach the independent particle limit Fq(x) — > as x — > 0. The simpler explana- 
tion is that the extrapolations to zero concentrations are too uncertain. For more accurate 
statements about the agreement (or disagreement) between theory and experiment, one 
should have many more susceptibility data especially at low concentrations. 

The next step in the theoretical procedure would be the summation of properly antisym- 
metrized exchange diagrams as suggested by Babu and Brown.0 But in view of the above 
discussion and the very high accuracy that is needed to extract the first antisymmetric 
Landau parameter from susceptibility data it might well be worthwhile to reconsider these 
experiments. 
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APPENDIX: DERIVATION OF THE SELF ENERGY 

In this Appendix we give the derivation of the self energy starting from the continuity 
equations ( |2.12| ) and ( |2.21| ). The first continuity equation defines the self energy in terms 
of the solution of the second continuity equation. We explain our method for solving these 
equations in momentum space leading to Eqs. ( |2.5| ) and ( |2.6| ) and use the notations of Ref. 



The impurity pair- and triplet distribution functions are defined in the usual way, 

^ (34 Hro,r 1 ) = -^ y p( 34 )(r ,r 1 ), (Al) 

S (344) (ro, ri ,r 2 ) = p^(r , ri ,r 2 ) , (A2) 
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and the impurity structure function S^(k) is the Fourier transform of the pair-distribution 
function g^ 34 \r) — 1. Analogous definitions are used for the background 4 He distribution 
and structure functions. 

Since our background 4 He liquid is homogeneous the continuity equations are easiest to 
solve in Fourier space. We have defined the Fourier transform of the fluctuating one-particle 
density in Eq. ( [2.11 ), for the time-dependent two-particle correlation function it has the 
following form: 



«gS(p) 



p (3) p (4) 



J d 3 r d 3 n dt e -i(kTo+P'(r - ri )-^) ^(34)^ ^.j) _ 



Using these notations the first continuity equation ( |2.5| ) transforms into the form 



(A3) 



2m?! 



<5p( 3 )(fc, 



(A4) 



with the self energy 

£ (3) (fc,w) 



d 3 p 



2m 3 J (27r) 3 p( 4 ) 



a 



(34), 



(A5) 



The second continuity equation (|2.6| ) contains the triplet distribution function. In order 
to calculate its Fourier transform we define the triplet structure function 



£(344)/ 



r ,ri,r 2 ] 



^ 344 )(ro,r 1; r 2 ) 



9 {34) (r ,r 2 ) 



9 



(44), 



n,r 2 ) + 2. (A6) 

Our approximation for this quantityi'i includes the convolution approximation together with 
the triplet correlation function. In the momentum space that leads to the result 

£( 344 )(k 1; k 2 , k 3 ) = 5(ki + k 2 + k 3 ){^ 34 )(k 2 ) SW(k 3 ) + 5( 34 )(k0 [^( 44 )(k 2 ) S^\k 3 ) - 1 
+ 4 344) (ki, k 2 , k 3 ) [S^XkJ + l] S^(k 2 ) S( 44 )(k 3 )} . (A7) 



The second continuity equation (|2.6| ) transforms in the momentum space into the following 
integral equation 



[hu - t^(k + p) - e^\p)} a£J(p) = hu 



where the kernel Kk }U1 (p, q) is 



k p S^jp) 
k 2 S(^(p) ' J (2tt)V 4 ) 



d 3 q 



"k 3 ?(q) ^(p,q) 



(A8) 



^(P,q) = S^\q){[(S^\\p - q|) + l) uf M \p - q, -p, q) + ^ 34 )(| P - q|) 



(k + p) ■ (k + q) 

p2 



t (3) (p) 



}-^ 34 )(| P -q|)^ e ( 4 )(p). (A9) 



The singularity structure of the self energy as well as the integral equation is best ac- 
counted by introducing the following notation, 
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f£$(p) = I** ~ i(3) ( k + P) - e(4) W] 4 3 S(P) • 
Inserting this into Eqs. ( |A5| ) and ( ]AS| ) we get our final forms for the self energy 

h 2 r d* P k'j>sM(p)p£2{j>) 



UJ 



2 m 3 7 (27r) 3 p {4) [^-t( 3 )(k + p) -e( 4 )(p) 



(A10) 



(All) 



and the integral equation 



d 3 q 



fc 2 S«*)(p)~ 7 (27T)3p(4) f^_t(3)( k + q )_ e (4) (?) 



(A12) 



The angle integration can be done conveniently by expanding /3^(p) in terms of Leg- 
endre polynomials, 



/Op)=EA (34J (^^K 



(34). 



(A13) 



1=0 



Here CC fep — COS 9k p with the angle 6k p between vectors k and p. Similarly we expand the 
other angle dependent quantities 

oo 

S( 34 )(|p-q|)^]>>^,g)P,(^), 



p q 

pq 

„~,( 344 ) 



1=0 

oo 



S^(\p- q \) = ^i(p,q)PiM 



p-q|) + i)4 (p-q>-p>q) =J2 s 'i(P: q) p i( x P q) > 



p • q 

pq 



(S (34) (|P - q|) + l) 4 344) (P " q, -P, q) = E *i(P. ?W*w) 



(A14) 



z=o 



Using these notations we can perform the angle integration analytically in terms of 
Legendre polynomials and the complex Leg endre functions of the second kind Qi(z), 



Qi(z) 



1 r \ x P ^ 



(A15) 



2 J -i z — x 

That leads to the coupled integral equation for the expansion coefficients flj (k, u,p) 



Pl u \k,u,p) + 



1 



S^^ihus^q) 



e{p) 



k Si(p,q) -pqai(p,q)} - -—qai(p,q) 



(A16) 
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Energies are expressed in units which include the effective mass, 

2m* H huj 



2 



fW = 2^ (ai7) 

and we have also introduced notations 

hou — k 2 — q 2 — e(q) 
Zkq= 2kq 
s t (p,q) = si(p,q) + s'^q) 

ai(p, q) = ai{p, q) + a[(p, q) . (A18) 
The functions T\ m (z kq \ Xf m (z kq ) and Tf m (z kq ,p, q) are the results of the angle integration 

T i /_ n _ / Pm(z kq )Qi(zk q ), when m < I , 
MrnVk*) | Q m {z kq ) Pl {z kq ), when m > I ' [ ^ } 

( z kq P m (z kq )Qi(z kq ), when m < I 
^irrX z k q ) = S z kq Q m (z kq )Pi(z kq ), when m > I , (A20) 

[ z kq Pi(z kq )Qi(z kq ) - 1/(21 + 1), when m = I 

\ ( I 1 n\ 2 _ , s. j P m (z kq )Q n (z kq ), when m < n , . 



with the 3j-symbol 



Z 1 n 
x J' 

Using the above notations the self energy has the form 
m* H f°° p 2 dp 



^ 3 \k,u) = ^ £^S^\p)Y,m",P)[z kp Qi(z kp ) -S^] . (A22) 

The analytic properties of the self energy as well as the integral equation are buried now 
into the Qi(z) function. It has a logarithmic singularity at \z\ — 1, it is a real function 
when 1 2 1 > 1 and complex when \z\ < 1. The singularity is integrable and thus the q-space 
integration in Eq. ( |A16|) and the p-space integration in Eq. ( |A22|) can be performed. 
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TABLES 



TABLE I. Pressure dependence of the hydro-dynamic effective mass from various calculations 
and experiments. The second column contains the result of our microscopic calculation described 
in sections [l| and IV A ; the next column contains the hydrodynamic effective mass as obtained 
from the fit (|4.8|) proposed by theory to the experiments of Ref. 28. Cols. 4 and 5 contain the 



hydrodynamic effective mass as obtained from the linear extrapolation and from the fit (4.8) to 



the experiments of Ref. 29 and Col. 6 the results from the fit to the magnetic susceptibility data. 



P (atm) 


This work 


Ref. ^8| (extrap.) 


m* H /m 
Ref. H 


Ref. ^ (extrap.) 


magn. susc. 





2.09 


2.18 


2.23±0.02 


2.15 


2.27 


5 


2.22 


2.31 








10 


2.34 


2.44 


2.52±0.02 


2.39 


2.42 


15 


2.45 


2.54 








20 


2.55 


2.64 


2.70±0.03 


2.62 


2.58 
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TABLE II. The half width of the peak in the dynamic structure function of the impurity 
excitation in the roton region for 1.8A -1 < k < 2.4A -1 . The third column gives the distance the 
impurity excitation can propagate within its life time. 



k A" 1 half width (K) distance (A) 

1.8 028 35J 

1.9 0.63 16.2 

2.0 1.01 10.3 

2.1 1.55 6.84 

2.2 2.3 4.67 

2.3 3.1 3.51 

2.4 4.0 2.74 
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TABLE III. Pressure dependence of the coefficients of the expansion ( |4,8| ) for the concentration 
dependence of the effective mass. The expansion coefficients a, b, c, and d are from Ref. 51, 



P (atm) abed 

L49 L39 -18.2 36/7 

5 1.07 3.00 -22.6 40.2 

10 0.789 4.48 -28.2 50.4 

15 0.501 6.17 -36.1 66.8 

20 0.310 7.41 -42.1 80.1 
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TABLE IV. Pressure dependence of the parameters of the fit ( 4.10 ) of the un-normalized Fermi 
liquid parameter {m/m*)FQ as obtained from the dynamic calculation. 



P (atm) abed 

0.447 -4.371 14.67 -22.35 

5 0.394 -3.710 10.82 -14.73 

10 0.362 -3.319 8.410 -9.660 

15 0.344 -3.012 6.336 -4.927 

20 0.326 -2.733 4.383 -0.349 
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FIGURES 




FIG. 1. The figure shows the theoretical hydro dynamic mass from our calculation described in 
Section || (solid line), and our zero-concentration extrapolations ( |4.6| ) and ( [4.7] ) of the data of Ref. 
p8| (long dashed line) and p9| (short dashed line). 
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FIG. 2. The momentum dependent part of the hydrodynamic effective mass for densities be- 
tween 0.022 A and 0.025 A as labeled in the fi gure. To highlight the momentum-dependent 
part, we display, according to Eq. (}4.2j), the ratio (m* H (k) — m* H (0))/(k 2 m* H (0)). 
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FIG. 3. The excitation spectrum of the 3 He impurity. The solid curve is the result of the 
present theory. It is compared with the measurements by Greywall (Ref. |3^) (long dashed line), 
Fak et al. (Ref. |30|) ( short dashed line) and Owers-Bradley et al. (Ref. |3^) (dash-doted line) . The 
dotted line shows, for reference, the experimental phonon-roton spectrum.E3 
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FIG. 4. The pole strength of the elementary impurity excitation mode is plotted as a function of 
momentum (solid line). The measured strength of the particle- hole excitation at 1% concentration 
and saturation vapor pressure from Fig. 11 of Ref. |3^ is shown with circles. For comparison we 
also show the effective mass as a function of momentum. 
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FIG. 5. The impurity dynamic structure function S^ 3 \k,uj) plotted in the k,co plane. Also 
shown are the phonon-roton spectrum of the background 4 He (heavy solid line, the data are from 
52), and the decay threshold of the impurity excitation mode (dashed line). 
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FIG. 6. The dynamic structure function of the impurity S^(k,u) in the roton region as a 
function of u; for fe=1.9A _1 , 2.1 A" 1 and 2.3 A" 1 . 
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FIG. 7. Theoretical and experimental effective mass ratio m as a function of pressure P 
and concentration x. The full curve is the fully self-consistent result, the dotted curve is the result 
without retardation effects, i.e. using the same quasiparticle interaction as in the spin-channel. 
The short dashed curve the static approximation. Symbols with error bars refer to the data of Ref. 
28. 



39 




40 




FIG. 9. The figure shows a comparison between experimental (diamonds, from Fig. 1 of Ref. 
39) and theoretical magnetic susceptibilities. The vertical scale is arbitrary and taken from Ref. 
39. The theoretical results were scaled to generate the best overall fit to the data at 0.27% and 
1.33% concentration. The solid lines are the results of the full microscopic calculation, the long 
dashed lines the results from using the static interaction W e g(q). Note that the static results for 
5.1 % concentration are close to the experiments at saturation concentration, and the static results 
for 8.8 % percent are off-scale. 
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FIG. 10. The figure shows the calculated Landau parameter Fq determined from the micro- 
scopic interactions V c s(q,0) (solid lines) and W e s(q) (dashed lines) together with the results of 
Ahonen et al. (boxes). The highest theoretical curves correspond to zero pressure and the lowest 
to 20 atm. Also shown are estimates for the Landau parameter Fq from the data of Ref. |39|, using 
the effective masses from Ref. |28| (diamonds). The lowest set of symbols corresponds to saturation 
pressure, the medium set to p = 10 atm, and the upper set to p = 20 atm. There is no experimental 
value at zero pressure and 8.8% concentration. 
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FIG. 11. The figure shows the coordinate-space effective interactions V e ff(r,0) (solid lines) 
and Weff(r) (dashed lines) for densities p = 0.020 A below saturation, p = 0.022 A close 

A_3 
close to solidification. The potentials with the higher 

repulsive peaks and deeper attractive wells correspond to the higher pressure. The Aziz-potentialll 

is shown for reference (dotted line). 
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FIG. 12. The figure shows the coordinate-space effective interactions V sca t(r) entering the 

0.022 A" 3 (solid line), p = 0.024 A" 3 (long-dashed 



scattering equation ( 4.12| ) at the density p 
line) and p 
line) . 



0.026 A (short-dashed line). The Aziz-potentialU is shown for reference (dotted 
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FIG. 13. The 3 He- 3 He scattering phase-shifts are shown at 1 percent concentration for the 
densities p = 0.022 A" 3 (solid line), p = 0.024 A -3 (long-dashed line), and p = 0.026 A~ 3 
(short-dashed line) for I = 0, 1, and 2. 
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FIG. 14. The 3 He- 3 He scattering amplitudes corresponding to the phase shifts shown in Fig. 
(|l~3|) at 1 percent concentration for the densities p = 0.022 A (solid line), p = 0.024 A 
(long-dashed line), and p = 0.026 A (short-dashed line) for I = 0, 1, and 2. The stars show 
the results in the simple HNC/0 approximation used by Owen@. 
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